First-principles and semi-empirical van der Waals study of thymine on Cu(llO) surface 
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In this study we investigated by means of density functional theory calculations the adsorption 
geometry and bonding mechanism of a single thymine (C5H6N2O2) molecule on Cu(llO) surface. In 
the most stable energetic configuration, the molecular plane is oriented perpendicular to substrate 
along the [110] direction. For this adsorption geometry, the thymine molecule interacts with the 
surface via a deprotonated nitrogen atom and two oxygen ones such that the bonding mechanism 
involves a strong hybridization between the highest occupied molecular orbitals (HOMOs) and the 
(i-states of the substrate. In the case of a parallel adsorption geometry, the long-range van der 
Waals interactions play an important role on both the molecule-surface geometry and adsorption 
energy. Their specific role was analyzed by means of a semi-empirical and the seamless methods. 
In particular, for a planar configuration, the inclusion of the dispersion effects dramatically changes 
the character of the adsorption process from physisorption to chemisorption. Finally, we predict the 
real-space topography of the molecule-surface interface by simulating scanning tunneling microscopy 
(STM) images. From these simulations we anticipate that only certain adsorption geometries can 
be imaged in STM experiments. 

PACS numbers: 71.15.Dx 68.43.-h 68. 43. Be 68.43.Fg 



I. INTRODUCTION 

The development of electronic devices like organic 
light-emitting diodes (OLEDs)ii^ has opened the road 
to an alternative technology to the silicon-based one to 
manufacture integrated electronic devices built on single 
molecules. Nowadays, the emerging field of molecular 
electronics^^^ focuses on the possibility of using organic 
molecules adsorbed on surfaces as basic functional units 
to construct technological relevant electronic devices 
such as single molecule diodes ^i^i^ organic field-effect 
transistors^A^ii^ or ultra- high-density memory circuits 

A very appealing feature of molecular devices is the 
prospect of tuning their physical properties at atomic 
scale by chemical functionalization of the molecule- 
surface interface under consideration. A key aspect of 
this approach is to employ molecules with several func- 
tional groups that adsorb on surface through one such 
group while the others can be used to react with other 
molecules to selectively modify the geometry and elec- 
tronic properties of adsorbate-substrates systems. For 
instance, experimental studies aimed to investigate the 
functionalization of the Cu(llO) surface by adsorption of 
organic molecules have been performed for terephthalio^ 
and oxalic^^ acids. 

In this context, the chemical functionalization of sur- 
faces via biologically relevant molecules is of particular 
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interest due to the potentially important applications 
such as molecular biosensors ^^i^i^ In particular, one 
of the most important molecules in biology is thymine 
(C5H6N2O2). Its importance stems from the fact that 
it is a fundamental nucleobase of the deoxyribonucleic 
acid (DNA) and encodes genetic information by pairing 
with the adenine base. The adsorption of thymine on 
surfaces is then of particular interest since it represents 
a possible channel for anchoring of DNA on substrates. 
The relevance of this process is also highlighted by the 
observation that the thymine base is absent in the case 
of the ribonucleic acid (RNA). 

From experimental point of view, the adsorption 
process of thymine on Cu(llO) surface was inten- 
sively investigated by means of reflection-adsorption 
infrared spectroscopy (RAIRS)r^i^ (soft-)X-ray pho- 
toelectron spectroscopy (XPS)r^i2£ near-edge X-Ray 
adsorption fine structure (NEXAFS)^i2£ and photo- 
electron diffraction^^ at room or higher temperatures. 
These studies suggest that, in the limit of low coverage, 
the geometry of the thymine-Cu(llO) system corresponds 
to an upright orientation of the molecular plane aligned 
along [110] surface direction. However, the chemical iden- 
tity of the functional groups involved in the bonding with 
the metal substrate is less clear. One possibility is to con- 
sider that the thymine molecule binds to Cu(llO) surface 
via a (dehydrogenated) N atom adjacent to two carbonyl 
groups (CO) while another possible adsorption geome- 
try implies a adsorbate-substrate interaction through the 
imino (CNH) and carbonyl (CO) groups. 

In consequence, in the present study we focus on an ab 
initio investigation of the adsorption geometry and bond- 
ing mechanism of a single thymine molecule on Cu(llO) 
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surface. The bonding geometry of the adsorbed molecule 
was analyzed as a function of the orientation of the molec- 
ular plane with respect to substrate. More specifically, 
we considered the case when the heterocyclic ring is per- 
pendicular to surface and aligned along its [001] and [110] 
directions. We also examined a parallel adsorption geom- 
etry such that the molecular plane lies flat on surface. In 
particular, for this bonding geometry the long-range van 
der Waals interactions are important from geometrical 
and energetic point of view. The effect of the disper- 
sion corrections on the adsorption geometry was taken 
into account by means of a semi-empirical approach while 
their impact on the adsorption energy was also evaluated 
within the seamless (vdW-DF) method. Besides this, the 
bonding mechanism of the thymine-Cu(llO) system in 
the perpendicular adsorption geometry was investigated 
considering a single and double deprotonated C5H6N2O2 
at N sites. In each case, we considered two scenarios 
for the bonding functional groups as already discussed 
above. The experimental fingerprints of all adsorption 
geometries analyzed in our study as revealed by scanning 
tunneling microscope (STM) was also simulated. 

II. THEORETICAL METHOD 

Our ab initio study was performed within the frame- 
work of density functional theory (DFT)^^ by using 
the generalized gradient approximation (GGA) for the 
exchange-correlation energy functional in the form pro- 
posed by Perdew-Burke-Ernzerhof (PBE).^^ The Kohn- 
Sham equation^ have been solved self-consistently us- 
ing the pseudopotential method24 as implemented in the 
VASP program j^^i^^ The electron-ion interactions are re- 
placed by pseudopotentials described by the projector 
augmented- wave method (PAW);^ 

The thymine-Cu(llO) system was modeled by a peri- 
odic slab geometry using the theoretical lattice param- 
eter of the bulk Cu (3.64 A). In each supercell, the slab 
consists of five atomic layers separated by a vacuum re- 
gion of^ 21 A. The size of the in-plane surface unit cell 
was set to 5x4. For this geometrical setup, the Kohn- 
Sham orbitals were expanded over a plane- wave basis set 
which includes all plane waves up to a cutoff energy Ecut 
of 450 eV. The geometry of the thymine-Cu(llO) sys- 
tem was optimized by relaxing the atomic positions of all 
molecule atoms and those in two surface layers. The equi- 
librium geometry of the molecule-Cu(llO) surface was 
obtained when the calculated Hellmann-Feynman forces 
were smaller than ^ 0.001 eV/A. Finally, the Brillouin 
zone integrations were carried out using only the F-point. 

To investigate the structural stability of different con- 
formational geometries of a single or double deprotonated 
adsorbate-surface systems, we calculated the adsorption 
energy Eads and adsorption enthalpy AH^^q at a tem- 
perature T of K for each geometry under consideration. 

The adsorption energy Eads is given by 

Eads = Egys — Ecu{110) ~ Eraolecule- (1) 



where Egys represents the total energy of the adsorbed 
thymine-Cu(llO) surface system, £^cn(iio) is energy of 
the isolated Cu(llO) surface and Emoiecuie denotes the 
energies of the isolated thymine molecular species as fol- 
lows: (i) EThymine -thymine and (ii) E^l^^.^^ -single 
(superscript H ) or —double (superscript H2) deproto- 
nated thymine molecule, respectively. 

The adsorption enthalpy AH^^q is defined as the dif- 
ference between the sum of the product energies {Egys 
and Eh2 —the total energy of the isolated hydrogen 
molecule H2) and the sum of reactant energies {Ecu{iio) 
and EThymine)- Therefore, the adsorption enthalpy 
AH^^Q is defined by 

AH^% ={Esys + factor * Eh,)- 

+ ^Cn(llO))- 

where factor is equal 0.0 (for thymine molecule), 0.5 
(for a single deprotonated case) or to 1.0 (for a double 
deprotonated thymine molecule). A negative adsorption 
enthalpy indicates an exothermic reaction between the 
thymine molecule and the Cu(llO) surface while a pos- 
itive value of AH^^Q implies that such reaction is en- 
dothermic and requires external energy to make it possi- 
ble. 

However, the double deprotonation process of the 
C5H6N2O2 molecule can take place through an inter- 
mediate step involving a single deprotonated thymine 
molecule. Two different scenarios are conceivable: (i) the 
single deprotonated molecule adsorbed on the Cu(llO) 
surface loses one of its proton such that the adsorption 
enthalpy is given by 

AH^% = «| + 0.5 * i?^^, ) - . (3) 

and (ii) a single deprotonated thymine species in gas 
phase adsorbs on the metal substrate where undergoes 
a second deprotonation such that 

AH^%={E^yl^0.5^EH,)- 

{^Thymine + ^Cn(llO))- 

Thus, the adsorption enthalpy AH^^q calculated for a 
double deprotonated C5H6N2O2 molecule by Eqs. [21 [3] 
and [Hwill allow us to get a hint on the details of the de- 
protonation process of this molecule on Cu(llO) surface. 

An important question arising when studying the ad- 
sorption process of molecules on metal surfaces concerns 
the role of long-range van der Waals interactions on the 
geometry and the adsorption energy of the adsorbate- 
substrate system in question. We addressed this prob- 
lem in our study by investigating the thymine molecule 
on Cu(llO) surface also in the framework of the semi- 
empirical approach (DFT-D) as proposed by Grimme.^^- 
The basic advantage of this method is the possibility to 
add the self-consistent calculated van der Waals forces to 
the Hellmann-Feynman ones to fully relax the molecule- 
surface system. We previously emphasized the impor- 
tance of the inclusion of these dispersion forces to obtain 
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a reliable adsorbate-substrate geometry in the case of the 
pyridine adsorbed on Cu(llO) and Ag(llO) surfaces^^ or 
in the case of a N,N/-di(n-butyl) quinacridone monolayer 
on the Ag(llO) substrate. A similar semi-empirical ap- 
proach was used to study the interaction of adenine 
with graphite(OOOl) surface^ or applied to noble-gas, 
N2 and benzene dimers.^^ Additionally, the energetics of 
the thymine-Cu(llO) substrate was analyzed by means 
of the seamless approach (vdW-DF) proposed by Dionet 
ali^ In this ab initio method, the vdW interactions are 
evaluated in a non-self-consistent procedure as a nonlocal 
correlation energy functional. This approach was already 
employed to investigate, for instance, benzene, naphtha- 
lene and phenol on graphite (000 1)^^^*^ or thiophene on 
Cu(llO) surfaces^ 

III. RESULTS AND DISCUSSION 
A. Thymine in gas phase 

To analyze how the molecule-substrate interac- 
tion modifies the geometry of the adsorbed thymine 
(C5H6N2O2) on Cu(llO) surface, we first optimized the 
geometry of this molecule in the gas phase [see Fig.[Ha)]. 
The structural optimization was performed by using a 
cubic unit ceh of 20x20x20A^. The calculated bond 
lengths and bond angles are reported in Table [H For 
comparison, in this table we also list the data of similar 
DFT calculations performed in Ref. [s^- Overall, one can 
observe that both studies provide a similar description of 
the molecular geometry. 

Similarly, the first step in understanding the electronic 
structure of the thymine-surface interface is to analyze 
the electronic structure of the C5H6N2O2 molecule in 
the gas phase. Therefore, in Fig.^a) we present the cor- 
responding angular-momentum resolved local density of 
states (LDOS). The highest occupied molecular orbital 
(HOMO) as well as the lowest unoccupied molecular or- 
bital (LUMO) are 7r-like molecular orbitals, originating 
mainly from the atomic p orbitals perpendicular on the 
molecular plane (pyrimidine ring, consisting of C and N 
atoms). The molecular orbital HOMO-1 and HOMO-2 
have a a character indicating that the chemical bonding 
takes place in the molecular plane. 

B. Ground-state geometry and energetics 

1. Perpendicular adsorption geometry 

Without vdW interactions - To explore the possibil- 
ity to chemically functionalize the Cu(llO) surface via 
the adsorption of a single thymine molecule on it, in our 
study we investigated several perpendicular adsorbate- 
substrate geometries. In all these adsorption configura- 
tions, the thymine binds to surface via the N atom after 
a deprotonation process involving a N-H bond breaking. 




(b) (c) 




(d) (e) 

FIG. 1: (a) Ball- and- stick model of the thymine molecule in 
the gas phase, (b) In the ground-state, the thymine molecule 
adsorbs on the Cu(llO) surface via three bonds (Nl-Cu,01- 
Cu,02-Cu) with the molecular plane perpendicular to sur- 
face and oriented along the [110] surface direction. This ad- 
sorption configuration is labeled as ONO-H in Table HT] An 
additional adsorption geometry obtained by removing the hy- 
drogen atom colored in grey (H2 in (a)) was also investigated 
(denoted as the geometry ONO-H2 in Table lUl. (c) Alter- 
natively, the thymine can attach to Cu(llO) surface via two 
chemical bonds (N2-Cu and 01-Cu). For this configuration 
one can consider also the cases of a single (NO-H geome- 
try in Table lU or a double (NO-H2 in Table |ll| deproto- 
nated thymine molecule, (d) Top and (e) lateral view of the 
most stable configuration of a parallel geometry of thymine 
on Cu(llO). The 01-02 "bond" line is oriented parallel to 
[001] surface direction. The plots presented in the present 
work have been obtained by using XCrysDen.— 



A binding of thymine on the Cu(llO) substrate through a 
N atom is similar to the pyridine on Cu(llO) and Ag(llO) 
surfaces. 

However, in the case of the thymine-Cu(llO) system, 
there are two different possible bonding N atoms denoted 
as Nl and N2 in Fig. [T](a). The molecule-surface geome- 
try corresponding to the adsorption through the Nl atom 
is presented in Fig.^b) (hereafter referred to as ONO-H 
configuration) and that corresponding to the bonding via 
the N2 atom is depicted in Fig. [T](c) (hereafter referred 
to as NO-H configuration). In both cases, the molecu- 
lar plane given by the pyrimidine ring is oriented along 
the [110] surface direction. In the present work we also 
investigated the case when the molecular plane is paral- 
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Gas phase Cu(llO) 
this work Ref. 37 



Configuration Perpendicular Parallel 

Orientation [110] [110] [001] 

DFT DFT+vdW 







Bond Length 


(A) 






Cl-01 


1.232 


1.222 


1.266 


1.232 


1.229 


C2-02 


1.229 


1.218 


1.262 


1.248 


1.251 


Cl-Nl 


1.412 


1.408 


1.376 


1.427 


1.430 


C2-N1 


1.385 


1.386 


1.356 


1.372 


1.369 


C2-N2 


1.392 


1.390 


1.379 


1.375 


1.370 


C3-N2 


1.375 


1.380 


1.367 


1.378 


1.379 


C3-C4 


1.358 


1.352 


1.360 


1.363 


1.362 


C4-C5 


1.511 


1.501 


1.511 


1.508 


1.506 


C4-C1 


1.479 


1.468 


1.468 


1.473 


1.475 


N-Cu 






1.988 






01-Cu 






2.017 






02-Cu 




Bond Angle 


2.028 


2.218 


2.135 


01-Cl-Nl 


119.9 


120.4 


120.8 


119.6 


119.6 


02-C2-N1 


124.5 




124.8 


123.8 


122.9 


C1-N1-C2 


128.6 


128.3 


121.6 


127.8 


127.6 


C2-N2-C3 


124.2 


123.9 


122.8 


123.1 


123.2 


N2-C3-C4 


122.9 


122.7 


121.1 


123.2 


122.8 


C3-C4-C5 


123.8 


124.1 


123.1 


123.5 


123.8 


C3-C4-C1 


117.8 


118.2 


116.7 


117.7 


118.0 


C4-C1-01 


125.9 




119.5 


126.5 


126.5 


C4-C1-N1 


114.2 


114.5 


119.7 


113.9 


113.8 



TABLE L The bond lengths and bond angles calculated for the thymine molecule in gas phase and adsorbed on Cu(llO) surface. 
In the case of a perpendicular adsorption geometry, only the values obtained for the ground state [see Fig. [T]^b)] are presented. 
In the case of a parallel adsorption geometry, only the data of the most stable configurations are reported. In particular, the 
geometry corresponding to the orientation [001] is sketched in Fig. [T]^d). In the orientation [110], the 01-02 "bond" line is 
oriented parallel to this surface direction. The atomic labels are those indicated in Fig. [TJa). 



lei to the [001] surface direction. Additionally, for each 
adsorption geometry already mentioned, we considered 
the possibility that the C5H6N2O2 molecule on Cu(llO) 
surface is double deprotonated, i.e., both N-H bonds are 
broken. 

To determine which adsorption configuration is ener- 
getically the most favorable one, we compared the corre- 
sponding adsorption energies reported in Table [III One 
can easily observe that, for all perpendicular thymine- 
surface configurations considered in our study, in the 
ground-state adsorption geometry the molecular plane is 
oriented parallel to [110] direction. Also, a lower adsorp- 
tion energies is obtained for the ONO-H adsorption con- 
figuration than for the NO-H one and a similar conclu- 
sion is drawn for the ONO-H2 geometry when compared 
with the NO-H2 one. This behavior is not surprising 
since early studies of formate^^ and terephthalic acid^^ 
on Cu(llO) surface revealed that the oxygen atoms of 
these molecules adsorb on top of the surface atoms. As 



clearly shown in Fig. [IJb) and (c), in the case of the 
[110] adsorption geometries, both thymine oxygen atoms 
lie above the surface Cu ones. One the other hand, in 
the case of the ONO-H (ONO-H2) configurations there 
are two oxygen atoms above the surface one, in contrast 
to the case of the NO-H (NO-H2) geometries where only 
one oxygen atom interacts directly with the surface. 

Also, the calculated adsorption energies for the geome- 
tries with a double deprotonated thymine on the Cu(llO) 
substrate are lower than those obtained for the single de- 
protonated adsorption configurations. However, since for 
these adsorbate-surface geometries the number of atoms 
is different, it is not possible to compare directly the cor- 
responding adsorption energies to infer which configura- 
tion is the most stable one. However, by comparing the 
calculated adsorption enthalpies AH^^q one can deter- 
mine if the deprotonation process of thymine molecule is 
favorable from energetic point of view or not. In particu- 
lar, this information is of peculiar importance for a dou- 
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adsorption energy adsorption enthalpy 

Eads (eV) AH^% (eV) 

geometry orientation DFT(PBE) DFT-D vdW-DF DFT(PBE) DFT-D 

PBE revPBE 

Eg. [2] Eg. [3] Eg. g] Eg. [2] Eg. [3] Eg. g] 

Perpendicular 



ONO-H 


[110] 


-3, 


,666 


-4.011 


-4.078 


-3.613 


-0, 


,618 






-0.621 








[001] 


-3, 


,027 


-3.396 


-3.566 


-3.157 


+0, 


,022 






-0.006 






ONO-H2 


[110] 


-4, 


,643 


-4.994 


-5.091 


-4.623 


+0, 


,832 


+1.450 


-2.217 


+0.822 


+1.443 


-2.510 




[001] 


-3, 


,902 


-4.282 


-4.438 


-3.994 


+1. 


,573 


+1.551 


-1.475 


+1.535 


+1.540 


-1.800 


NO-H 


[110] 


-2, 


,674 


-3.050 


-3.042 


-2.668 


-0, 


,064 






-0.098 








[001] 


-2, 


,561 


-2.960 


-2.987 


-2.610 


+0, 


,048 






-0.007 






NO-H2 


[110] 


-3, 


,160 


-3.590 


-3.570 


-3.198 


+2, 


,316 


+2.380 


-0.294 


+2.228 


+2.326 


-0.578 




[001] 


-3, 


,136 


-3.510 


-3.591 


-3.208 


+2, 


,340 


+2.292 


-0.270 


+2.307 


+2.314 


-0.450 














Parallel 
















flat 


[110] 


+0, 


,296 


-0.260 


-0.516 


-0.002 
















[001] 


+0, 


,112 


-0.319 


-0.703 


-0.247 

















TABLE II: The adsorption energies Eads and adsorption enthalpies AH^^q for the thymine-Cu(llO) adsorption geometries 
considered in our study [see Fig. [TJb), (c), (d) and (e), respectively]. Here DFT(PBE) denotes the first-principles DFT results 
obtained with the PBE^^ exchange-correlation energy functional, DFT-D represents the DFT results plus the semiempirical 
vdW contribution as proposed by Grimme^ and vdW-DF implies the DFT results obtained by using the Dion^^ functional. 
Note that in the case of the vdW-DF calculations we employed PBB^ as well as revPBEj^ exchange- correlation energy 
functionals. 



ble deprotonated thymine since in this case the AHjcf^ 
and the relative surface free energies A7 plotted as a 
function of the chemical potential jj^h of the H atom will 
allow us to determine which reaction path corresponding 
to Eq. O Eq. [3]or to Eq. [4]is the most probable one. 

In the light of the above considerations, from Ta- 
ble [ni one can deduce that for a single deprotonated 
thymine molecule on Cu(llO) surface (ONO-H and NO- 
H configurations) the ground-state [110] adsorption ge- 
ometry does not require external energy to deprotonate 
the C5H6N2O2 molecule, in contrast to the [001] config- 
uration. 

In the case of a double deprotonated molecule, the 
large adsorption energy which is a measure of the binding 
strength suggests that the double deprotonated species 
might coexists together with single deprotonated species 
on the metal surface. An important question is how a 
double deprotonated species would arise. In our study 
we considered three different scenarios for a double de- 
protonation process: 

(i) the reaction path suggested by Eq. [2] implies an instan- 
taneous double deprotonation process when the thymine 
adsorbs on the Cu(llO) surface. This could take place 
at high temperatures since the reaction is endoterm (re- 
quires external energy to take place); 

(ii) the reaction path described by Eq. [3] implies that the 
single deprotonated thymine adsorbed on the Cu(llO) 
loses instantaneously its second hydrogen yielding a dou- 
ble deprotonated species adsorbed on the surface. This 
process is also an endothermic one and requires a larger 
amount of energy to take place as compared to (i); 



(iii) the Eq. [4] assumes a reaction path in which a sin- 
gle deprotonated thymine molecule in gas phase adsorbs 
on the metal substrate where it undergoes a second de- 
protonation process. Since a large amount of energy is 
released, the adsorption process is an exotherm one and 
the probability to occur is can be significant (depending 
on the hight of the energy barrier of the transition state). 

Since the adsorbed species of the flat thymine, ONO-H 
and ONO-H2 differ only through the number of hydrogen 
atoms, in Fig. [2] we have plotted the relative surface free 
energies A7 as a function of the chemical potential jj^h of 
the H atom. Details of how A7 was calculated are pre- 
sented in Appendix. At a given temperature, the chemi- 
cal potential of the H atom can be expressed in terms of 
the H2 molecular pressure which together with the tem- 
perature represent the parameters used in experiments to 
control which conformation of the thymine (single or dou- 
ble dehydrogenated) is adsorbed on the Cu(llO) surface. 
The most stable configuration is the one that minimizes 
the relative surface free energy A7. In consequence, the 
Fig. [2] indicates that, independent on the temperature, 
the most stable configuration is ONO-H for the chemical 
potentials of the H atom between —0.9 and 0.0 eV. Note 
that these values correspond to the chemical potential in- 
terval in which the experiments can be performed. How- 
ever, the experiment^ii^ii^i2£ indicate that by preparing 
the sample at 300 K the ONO-H species are present on 
the surface while when increasing the temperature higher 
than 500 K the ONO-H2 species are predominat at an 
experimental pressure of 10~^^ atm. Based on the the 
experimental observations and as shown by our ab initio 
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H pressure (atm) 




-0.5 



FIG. 2: Relative surface free energies A7 of the flat, ONO-H 
and ONO-H2 configurations as function of the H chemical 
potential /j^h- The values of the A7 are normalized by the 
surface area A of the unit cell (see Appendix). At the top of 
the X— axis, the chemical potential of the H2 molecule /JH2 is 
given as a pressure scale at two fixed temperatures: 300 K 
and 500 K. 



calculations (see Fig. [2j) we conclude that at room tem- 
perature (~ 300 K) there is a high activation energy be- 
tween ONO-H and ONO-H2 adsorbed species and there- 
fore the second deprotonation process cannot take place. 
Increasing the temperature over 500 K this energetic bar- 
rier can be overcome such that the ONO-H2 species are 
formed*^ 

To conclude, our first-principles calculations suggest 
that a double deprotonated thymine molecule adsorbed 
on Cu(llO) surface can be formed only at high tempera- 
tures as follows: (i) the thymine adsorbs on the Cu(llO) 
surface undergoing an instantaneous double deprotona- 
tion process (see Eq. [2])or (ii) the single deprotonated 
thymine adsorbed on the Cu(llO) loses its second hy- 
drogen yielding a double deprotonated species adsorbed 
on the surface (process described by Eq. [3]). Note that 
even if the adsorption enthalpy calculated for the depro- 
tonation process described by Eq. [4] is negative and thus 
this process is exothermic, a single deprotonated pyridine 
molecule in gas phase requires to be desorbed from the 
Cu(llO) surface. In consequence, to initialize this pro- 
cess one needs to provide an amount of externel energy 
that equals at least the adsorption energy of the ONO-H 
or NO-H configurations and therefore this scenario is less 
probable than those described by Eqs. [2] and [3l 

Role of the vdW interactions - When including the 
van der Waals interactions using the DFT-D approach, 
all relaxed thymine-Cu(llO) geometries considered in our 
study remain practically the same as those obtained with- 
out including the effect of the dispersion corrections. 



This outcome of our ab initio calculations is not sur- 
prising since in all perpendicular adsorption configura- 
tions the thymine molecule form a chemical bond with 
the surface via a N atom. As already discussed by Thon- 
hauser et a/.,^^ for covalently bonded systems like a CO2 
molecule or the bulk silicon, the inclusion of the vdW in- 
teractions does not alter the geometry obtained by DFT 
calculations. 

However, the dispersion corrections will generally lower 
the adsorbtion energy of an adsorbate-substrate system. 
The importance of this effect for chemisorbed systems 
was emphasized in the study of the adsorption mecha- 
nism of benzene on Si(OOl) surface performed by John- 
ston et al^ Their first-principles calculations based on 
the vdW-DF approacf>2^ revealed that the PBE and re- 
vised PBE (revPBE)^ calculations predict a different 
ground-state geometry than the vdW-DF ones. There- 
fore, to take into account the weak vdW interactions 
might be important even for the strongly chemisorbed 
molecule-surface systems. As presented in Table HH in the 
case of thymine-Cu(llO) surface the [110] adsorption ge- 
ometry remains the ground-state one for all four adsorp- 
tion configurations even in the presence of the attractive 
van der Waals interactions. In the semi-empirical (DFT- 
D) approach, their main effect is to lower the adsorption 
energy calculated for the perpendicular adsorbate-surface 
geometries by about 0.3 eV with respect to the values ob- 
tained with PBE. Note that this energy gain due to the 
dispersion corrections is similar to that obtained for pyri- 
dine on Cu(llO) and Ag(llO) surfaces. As regarding the 
adsorption enthalpies, their values are also lowered when 
the dispersion effects are included. 

A final note of this section concerns the comparison of 
the adsorption energies obtained with the semi-empirical 
(DFT-D) and seamless (vdW-DF) methods, the latter 
values being obtained with a code developed by us 4^ We 
briefly remind that in the vdW-DF approach the corre- 
lation part of the GGA functional is replaced by a sum 
of the LDA correlation functional Elda,c and a nonlocal 
term E^^: 



EvdW-DF = Egga - Egga,c + Elda,c + E\ 



ml 



(5) 



where Egga is the self-consistent GGA total energy 
Egga evaluated in a conventional DFT calculation. The 
nonlocal correlation energy E^^ can be expressed as a 
function of the charge density n(r) 



Ef = \ j j drdr'n{r)(^{r,r')n{r'). 



(6) 



where the kernel function 0(r, r ') is discussed in detail in 
Ref. 113. Note that in our study we employed both PBE 
and revPBE flavours of the GGA exchange-correlation 
energy functional. 

From Table [III it becomes apparent that in the case 
of the perpendicular adsorption geometries, the vdW- 
DF values calculated with PBE are in excellent agree- 
ment with those obtained by using the semi-empirical 



approach DFT-D. This observation leads us to the pre- 
Uminary conclusion that, in the case of the covalently 
bonded systems, the DFT-D and vdW-DF will provide 
similar results based on the PBE functional. However, 
when using the revPBE, the adsorption energies calcu- 
lated with the vdW-DF approach are higher by ^0.3 eV 
that those obtained with DFT-D method based on PBE. 
This behaviour is similar to that observed for chemisorp- 
tion energy of atoms and molecules on transition-metal 
surfaces as pointed out by Hammer et al^ 



2. Parallel adsorption geometry 

Role of the vdW interactions - In our study we also 
focused on the thymine-Cu(llO) surface interaction when 
the molecule lies flat on the substrate. More specifically, 
we considered two different adsorption configurations 
such that the molecule is oriented with the line given by 
0-N-O along [001] direction ([001] geometry) or along 
[110] one ([110] geometry). Intuitively, in these adsorp- 
tion configurations the molecule-surface bonding mecha- 
nism implies an interaction between a tt molecular orbital 
and surface states which eventually leads to a weakly 
bonded system. As was already discussed in literature for 
benzene and naphthalene molecules on graphite(OOOl)^ 
phenol on graphite(OOOl) and a-Al2O3(0001)^^ or for 
thiophene on Cu(llO) surfaces, the van der Waals cor- 
relations play a crucial role to correctly describe the en- 
ergetics of such systems. 

Indeed, the PBE calculations performed for both flat 
geometries predict an average molecule-substrate dis- 
tance of ^3 A A that is a clear fingerprint of a physisorp- 
tion process. However, by including the dispersion ef- 
fects and relaxing the adsorbate-substrate interface, in 
both configurations the molecule-surface distance drops 
to ^2.7 A. Besides this, the thymine becomes tilted with 
respect to the surface by an angle of ^111° and ^104° 
for [001] and [110] geometries, respectively. This behav- 
ior is similar to that observed for pyridine on Cu(llO) 
and Ag(llO) surfaces^^ and clearly emphasizes the im- 
portance of taking into account the structural relaxations 
for such weakly bonded molecule-surface systems. From 
bonding point of view, the main difference with respect 
to the perpendicular adsorption configurations is that a 
flat thymine molecule is anchored to Cu(llO) via a 0-Cu 
bond with a bondlength of ^2.1 A. 

The necessity to include the dispersion corrections 
when investigating van der Waals systems is illustrated 
by the analysis of the calculated adsorption energies. As 
shown in Table [TTl the PBE calculations predict actu- 
ally a non-bonding thymine-surface configuration since 
the adsorption energies are positive. However, when con- 
sidering the effect of the van der Waals interactions on the 
semi-empirical level, these energies are lowered on aver- 
age by ~0.5 eV and become negative, denoting that a flat 
adsorption geometry of thymine on Cu(llO) substrate is 
possible. In this case the [001] adsorption geometry is 




FIG. 3: (Color online) Perspective 3D-view of the the non- 
local correlation energy difference AE^^ of the thymine- 
Cu(llO) surface system with respect to its isolated surface 
and molecule components. The red color shows regions with 
a large contribution to the binding energy while the green 
color represents regions with a much smaller contribution. 

more stable than the [110] one. 

The use of the vdW-DF method does not qualitatively 
change the physical picture highlighted above. From 
quantitative point of view, the vdW-DF adsorption en- 
ergies obtained with the PBE functional are lower by 
a factor of ^2 with respect to those evaluated with the 
DFT-D approach. Similarly to the case of the perpendic- 
ular thymine-surface geometries, the vdW-DF adsorption 
energies calculated with the revPBE are about ^0.5 eV 
higher then their PBE counterparts. Although the rela- 
tive stability of the flat adsorption geometries considered 
in our study is not changed with respect to the PBE or 
the semi-empirical method, the use of revPBE leads to a 
very small value of the adsorption energy in the case of 
the [110] configuration. This is consistent with the gen- 
eral tendency of the revPBE functional to provide theo- 
retical chemisorption energies higher than the PBE one 
and thus presumably closer to the corresponding exper- 
imental values^ However, note that the revPBE does 
not generally lead to an overall better physical descrip- 
tion than the PBE functional. For instance, in the case 
of a typical van der Waals system such as Kr dimer— 
while the depth of the Kr-Kr interaction potential is bet- 
ter described by revPBE, the equilibrium distance of this 
dimer is better reproduced by the PBE functional. 

Since we applied the vdW-DF method to an initially 
weakly bonded system that in the final relaxed configu- 
ration exhibits a chemical bond, it is interesting to de- 
termine the role of the nonlocal correlation effects in this 
case. Therefore in Fig. [3] we present a plot of the dif- 
ference in the nonlocal correlation energy due to the ad- 
sorption process, i.e., we map in real space the quantity 

AEf =^n^s2ys_£:nZ,Cn(110)_^n/,t/iymme^ The peculiar 



8 




+1.9 eV 



-1.9 eV 



Energy (eV) 



Energy (eV) 



(e) 



(f) 



FIG. 4: Local density of states (LDOS) calculated for (a) 
the thymine in the gas phase, for (b) a flat thymine-surface 
adsorption geometry [see Figs. [TJd) and (e)]and for the (c) 
ONO-H, (d) ONO-H2, (e) NO-H and (f) NO-H2 perpen- 
dicular adsorption geometries of the C5H6N2O2 molecule on 
Cu(llO) surface. 




(i) 



(j) 



feature displayed by this image is a strong contribution 
of the nonlocal effects to the bonding energy arising from 
the charge density n(r) close to the O and Cu atoms that 
form a chemical bond. Qualitatively this behavior can be 
understood by noticing that to form a chemical bond be- 
tween the O and Cu atoms, the charge density will shift 
from these atoms to the space between them. Therefore, 
they become polarized which favors a significant contri- 
bution of the nonlocal correlation energy given by Eq. [6] 
to the binding energy of the thymine- Cu( 110) system. 



C. Electronic structure 

The calculated atom-projected local density of states 
(LDOS) for all adsorption configurations are presented 
in Fig. m A general feature of the electronic structure 
of all perpendicular adsorption geometries considered in 
our study is a strong hybridization of the highest occu- 
pied molecular orbitals (HOMOs) with the (i-states of 
the Cu(llO) surface in an energy range of 5 eV below 



FIG. 5: The simulated STM topographies for (a,b) ONO-H, 
(c,d) ONO-H2, (e,f) NO-H and (g,h) NO-H2 perpendic- 
ular adsorption configurations. The similar images for a flat 
adsorption geometry [see Fig. [U^d) and (e)] are shown in (i) 
and (j). The applied bias voltages are +1.9 and -1.9 eV, re- 
spectively. At a positive applied bias voltage the electrons 
will tunnel from the STM tip to the unoccupied electronic 
states of the molecule-surface interface. 



the Fermi level. This strong molecule-surface interaction 
leads to broad bands with a mixed molecular-metal char- 
acter. Besides this, in the energy range from -5 to -10 
eV, the interaction of the molecular states with s- and p- 
states of the substrate become also important. As regard- 
ing the most stable flat adsorption geometry depicted in 
Fig. [TJd), its electronic structure also involves mainly a 
strong hybridization of HOMOs with the d-states of the 
substrate. 

However, in an energy range of ±2 eV around Fermi 
level, the LDOS has specific features for each studied 
adsorption configuration. This observation is important 
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since the electronic structure of adsorbate-metal surface 
in this energy range can be probed by using the scanning 
tunneling microscope (STM). Therefore, using the Ter- 
soff and Hamann theory^ we simulated maps of constant 
current for the ground-state perpendicular and parallel 
adsorption geometries considered in our study. The the- 
oretical STM topographies are depicted in Fig. [5] for two 
different applied bias voltages. In general, the simulated 
STM images exhibit peculiar features for each adsorbate- 
surface configuration at each value of the bias voltage ex- 
cept for the NO-H and NO-H2 geometries which exhibit 
similar features. 

Among the simulated theoretical STM images, those 
obtained for the ONO-H geometry at a bias voltage of 
±1.9 eV present the most significant topographic dif- 
ference. Also, except for this configuration, the sim- 
ulated images suggest the possibility to obtain a sub- 
molecular resolution for single or double deprotonated 
thymine molecule adsorbed in a perpendicular configura- 
tion on Cu(llO) surface. For instance, the methyl group 
(CH3) appears as a distinct clear feature, especially in 
the case of the ONO-H2 and NO-H geometries although 
its geometrical form is specific for each configuration in 
question. A similar sub-molecular resolution of the STM 
images is observed for a parallel adsorption geometry cor- 
responding to the configuration shown in Fig. [TJd). 



IV. SUMMARY 

By using the ab initio pseudopotential method, we 
analyzed the adsorption process of a single thymine 
molecule on the Cu(llO) surface. To determine the 
ground-state adsorption geometry, we investigated sev- 
eral adsorbate-substrate configurations such that the 
molecular plane is perpendicular or parallel to the (110) 
surface. To explore the possibility to chemically func- 
tionalize this molecule-surface system, for the perpen- 
dicular adsorption geometry we considered the case of 
a single and double deprotonated thymine molecule an- 
chored to surface via different functional groups. Our ab 
initio simulations indicate that the strongest thymine- 
surface interaction takes place when the heterocycle ring 
is perpendicular to substrate and aligned along its [110] 
direction. For this ground-state adsorption geometry, the 
bonding of thymine on Cu(llO) surface is due to a strong 
hybridization of the highest occupied molecular orbitals 
(HOMOs) with the d-states of the substrate. In the case 
of a fiat adsorption geometry, our calculations demon- 
strate that the long-range van der Waals interactions 
play a very important role from both geometrical and 
energetic point of views. The effect of the dispersion cor- 
rections on the geometry of the relaxed adsorbate-surface 
interface was determined by using a semi-empirical ap- 
proach. If the relaxed geometry obtained from PBE cal- 
culations remains almost parallel to surface, the inclu- 
sion of the van der Waals interactions leads to tilted 
adsorption geometry by an angle of 104°. Moreover, 



the van der Waals corrections have a more dramatic im- 
pact on the calculated adsorption energies. The PBE 
adsorption energies obtained for the fiat adsorption ge- 
ometries are positive and thus the PBE calculations sug- 
gest a non-bonding molecule-surface configuration. On 
the contrary, the adsorption energies calculated with the 
semi-empirical and the vdW-DF methods are negative 
and thus the parallel adsorption configuration becomes 
bonding. In this case the inclusion of the van der Waals 
interactions significantly modifies the character of the ad- 
sorption process from physisorption to chemisorption. As 
a final remark, we also simulated the scanning tunneling 
microscopy (STM) topography of the thymine-Cu(llO) 
surface that can serve to distinguish in experiments be- 
tween various adsorption configurations of the thymine- 
Cu(llO) interface. 
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Appendix 

A detailed description of how to evaluate the surface 
free energies using the total energies obtained from DFT 
calculations can be found in Refs. 48 and i49|. In the 
following we will derive the equations used in the present 
work to calculate the relative surface free energies A7 
plotted in Fig. [2 

For a deprotonation process characterized by the rel- 
ative Gibbs free energy AG^^*, the relative surface free 
energy (A7) can be written as 

A7 = {1/A) * (AG^^" + factor * ^h)^ 

where A represents the surface area of the unit cell and 
jj.H2 is the chemical potential of the H2 molecule. The 
factor is equal 0.0 for thymine molecule, 1.0 for a single 
deprotonated or to 2.0 for a double deprotonated thymine 
molecule. 

Assuming that the entropic term is small as compared 
to the values of the adsorption entalpies at T = K,^^ we 
approximate AG^^^ by AH^^q. Therefore, the relative 
surface free energies A7 of the ground-state fiat, ONO-H 
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and ONO-H2 configurations are given by: where /i/f has been calculated according to Eq. (22) from 

the Ref. i49i and the values of AH^^q correspond to those 
A * A7^^^, ^ Ai^^i'o + 0-0 * Mi^ = -0.260 + 0.0 * /hh, obtained with the DFT-D method. 

A * A7^^^_^ ^ AH^% + 1.0 * = -0.872 + 1.0 * /i^, 
A * A7^^^_^^ ^ AH^Jq + 2.0 * = +0.544 + 2.0 * /i/,. 
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